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. A deep 75 ks Chandra ACIS-I data of NGC 2024 was analyzed, aiming at a 

search for diffuse X-ray emission in this one of the most nearby (415 pc) massive 
star-forming regions. After removing point sources, an extended emission was 
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detected in the central circular region with a radius of 0.5 pc. It is spatially 
associated with the young massive stellar cluster. Its X-ray spectrum exhibits 
a very hard continuum (kT > 8 keV) and a sign of He-like Fe K a line with the 
0.5-7 keV absorption corrected luminosity of 2xl0 31 ergs s _1 . Undetected faint 
^ | point sources, estimated from the luminosity function of the detected sources, 

contribute less than 10% to this emission. Hence the emission is truly diffuse 

O ! 

\q \ in nature. Because of the proximity of NGC 2024 and the long exposure, this 

discovery is one of the most strong supports for the existence of the diffuse X-ray 
emission in massive star-forming regions. 
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1. Introduction 

Increasing evidences of diffuse X-ray emission has been found in massive star-forming 
regions (MSFRs) with the Chandra X-ray Observatory; the Rosette Nebula (at a distance 
of D = 1.4 kpc; Townsley et al. 2003), M17 (D = 1.6 kpc; Townsley et al. 2003), RCW 
38 (D = 1.6 kpc; Wolk et al. 2002), NGC 6334 (D = 1.7 kpc; Ezoe et al. 2006), RCW 
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49 (D =1.9-7.9 kpc; Townsley et al. 2004), W51A (D =5.5-7.5 kpc; Townsley et al. 2004), 
NGC 3603 (D — 7 kpc; Moffat et al. 2002), the Arches Cluster (D = 8.5 kpc; Yusef-Zadeh 
et al. 2002), and the Quintuplet Cluster (D = 8.5 kpc; ; Wang et al. 2002). Townsley et 
al. (2003) explained the diffuse soft X-ray emission found in M17 in the context of strong 
shocks by fast stellar winds from young massive stars (Dyson and de Vries 1972; Castor et al. 
1975; Weaver et al. 1977). Recent results on NGC 6334 by Ezoe et al. (2006) indicate that 
the spectra of the diffuse emission varies from place to place; those in tenuous molecular 
cloud regions are soft and thermal with temperatures of several keV, while those in dense 
cores exhibit harder continua with a photon index of T ~1. They also have shown that 
these thermal and non-thermal spectra of the diffuse X-ray emission in MSFRs, found as a 
mixture in the NGC 6334 case, may be generally explained by the stellar-wind shock model. 
In spite of these observational progresses, even with Chandra, there remains an uncertainty 
how much undetected faint point sources contribute to the emission, because these MSFRs 
are relatively distant (> 1.4kpc). 

In order to unambiguously examine the diffuse X-ray emission in MSFRs as possible 
as we can, we analyzed archival Chandra data of NGC 2024. This region, known as the 
Flame nebula, is one of the nearest MSFRs (D = 415 pc; Anthony-Twarog 1982), located 
in the Orion Nebula. It is an HII region considered to be illuminated by an 08V-B2V star 
IRS 2b (Bik et al. 2003). In the vicinity of IRS 2b, there are one early B star candidate 
associated with an ultra compact HII region G206. 543-16. 347 and an infrared source IRS 
2 (Lenorzer et al. 2004), and seven compact dust condensations named FIR 1-7, possibly 
massive protostars (Mezger et al. 1988). In addition, ~ 300 low-mass ( ^ 2M ) young stars 
have been found by near infrared observations (Lada 1991). The estimated age of NGC 2024 
ranges from 0.3 (Meyer 1996) to several Myr (Comeron et al. 1996). A previous analysis of 
the same Chandra data on point sources has been published by Skinner et al. (2003). In this 
paper, we focus on a search for the diffuse X-ray emission in this representative MSFR. 



2. Observation 

Chandra observed NGC 2024 on 2001 August 8-9 using the ACIS-I for 21.9 hr. We 
started with the level 1 event files in the same way as Skinner et al. (2003). Utilized analysis 
software versions for the standard data reduction are different; we used the CIAO (Chandra 
Interactive Analysis of Observations) software version 2.3 and the calibration data base 
version 2.18. These new softwares allowed us to correct the data for the effect of the charge 
transfer inefficiency, while it was not possible in Skinner et al. (2003). No background flares 
were seen during the observation with the average count rate of the six ACIS chips of 9.1 ct 
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s 1 . We then excluded > 1.2 times of the average rate. This procedure excluded 2% of the 
exposure time. After these data screenings, the nominal exposure has become 75.3 ksec. 

We searched the data for diffuse X-ray emission, following to the analysis method of 
Ezoe et al. (2006). We first created adaptively smoothed X-ray images in two energy bands, 
as shown in figures 1 (a) and (b). We see a sign of an extended emission associated with 
the massive star IRS 2b and its vicinity. The extended emission in 0.5-2 keV is strong in 
the north-west direction of IRS 2b, while that in 2-7 keV is elongated in the north-west to 
the south-east direction. We then used the wavedetect program to identify sources using 
images in three energy bands (0.5-2, 2-8 and 0.5-8 keV) independently, in order not to miss 
very soft or strongly absorbed sources. The significance criterion and wavelet scales were 
set at 1 x 10~ 6 and 1-16 pixels in multiples of y/2. In the ACIS-I field of view, 301 sources 
were detected. Among them, 28 sources are newly identified ones, not listed in Skinner et al. 
(2003), due to our three-band searching method. 

We excluded all the detected point sources by creating a point source mask using the 
"Chandra Ray Tracer" 1 . For each source, we defined a radius to include ~ 98 % of photons 
at the Al Ka-line energy (1.497 keV). Then we excluded all these regions by applying the 
mask to the raw ACIS-I image, and created images of the residual emission using the CIAO 
tools dmf ilth and csmooth. Figures 1 (c) and (d) show thus obtained images of the extended 
emission. In order to evaluate the significance of this emission, we defined a circular region 
named CI (the large circle in the panel d of figure 1) with a radius of 4' =0.5 pc. The total 
area of the CI is 43/39 arcmin 2 or 0.62/0.52 pc 2 before/after excluding the area around 
point sources. The 0.5-7 keV count rate of the CI is 0.142 ± 0.001 cts s _1 , while that of the 
same region in blank-sky data is 0.084 ±0.001 cts s _1 . Errors are la statistical ones. Hence, 
its residual count rate is 0.057 ±0.001 cts s _1 or 4270±110 cts. It is thus clear that a highly 
significant extended emission is present in the central region of NGC 2024. 



3. Extended X-ray Emission 

In the presence of the significant excess emission, we immediately considered a possibility 
of the diffuse emission. To know its basic properties, we compared its background subtracted 
spectrum in figure 2 (a) with that summed over 176 point sources detected within the CI. 
The weighted ARF (ancillary response function) and RMF (response matrix function) were 
calculated using the CIAO programs mkwarf and mkwrmf , respectively. The apply_acisabs 
script was utilized when creating ARF files, to correct them for the decrease in the ACIS 
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quantum efficiency. The background spectra are extracted from the same regions in the 
blank-sky data for individual regions. From figure 2 (a), we can see important features of 
the extended emission. First, it is about one order of magnitude fainter than the summed 
point sources. Second, the extended emission show a harder continuum in the 1-7 keV range. 
Third, a hint of an emission line is seen in 6-7 keV. 

To know basic parameters of the extended emission, we then conducted spectral fitting. 
We employed a simple power-law model with an interstellar absorption and a narrow Gaus- 
sian. Here and hereafter all quoted errors in the spectral fitting refer to 90% confidence levels 
unless otherwise noted. The result is shown in figure 2 (b), where table 1 lists the obtained 
parameters. The fit was not acceptable with x 2 l u ~ 1-5 because of the excess around 0.5-1 
keV. The line center energy indicates a He-like Fe K a line from thermal plasma. Then, we 
consider an alternative "leaky absorber" condition; a single thermal emission component 
reaches us via two (or more) paths with different absorptions. This situation is possible in 
the MSFR like NGC 2024 where the density of the molecular cloud varies from place to 
place. We hence fitted the spectrum by a sum of two thermal components with independent 
absorptions, but with their temperatures and abundances together. The result is shown in 
figure 2 (c) and table 1. As a thin thermal plasma model, we utilized the APEC (astrophys- 
ical plasma emission code) 2 model. Then, the fitting result is improved from x 2 / u ~ 1-5 
to 1.3, which is significant with 99.7% confidence according to an F-test. The soft excess is 
represented by the mildly absorbed high temperature emission. The best fit temperature is 
high 11 keV, consistent with the small photon index of T ~ 0.9 in the power-law model fit. 
The 0.5-7 keV flux obtained from the leaky absorption model is 1.1 xlO -12 ergs s -1 cm s~ 2 . 

For comparison, we quantified the summed point-source spectrum. In the same manner 
of the extended emission, the source spectrum, the background, ARF, and RMF files were 
obtained. The fit for a single thermal emission model with one absorption was not acceptable 
with x 2 l v ~ 7 because this model cannot represent both soft excess below ~ 2 keV, similar 
to that in the extended emission case, and also various emission lines. Hence, we used the 
leaky absorption model with free Ne, Mg, Si, S, Ar, Ca, and Fe abundances, in order to 
represent the data better. We obtained results as shown in table 2 and figure 3. The best- 
fit temperature of 4.4 keV is consistent with typical values of young low-mass stars (e.g., 
Imanishi et al. 2001) and significantly lower than that of the extended emission. 

In spite of these spectral analysis, We here must evaluate an effect of photons escaping 
from the summed point sources because the point sources are far brighter than the extended 
emission. In the same way as Ezoe et al. (2006), we took into account the summed point 
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source spectrum by multiplying the best-fit model for the summed point sources by the 
energy-dependent escaping- fraction curve, estimated by the u Chandra Ray Tracer" (ChaRT). 
We have found that this contributes ~ 40% to the extended emission, we remove their 98% 
photons at 2 keV. After correcting the escaping photon effect, the 0.5-7 keV flux of the 
extended emission becomes 6.3xl0~ 13 ergs s _1 cm s -2 , yielding an absorption-uncorrected 
luminosity of 1.3xl0 31 ergs s _1 or the surface brightness of 2xl0 31 ergs s _1 pc~ 2 at an 
assumed distance of 415 pc. 

We refitted the CI spectrum including the escaping photon effect. The result is shown 
in figure 2 (d) and table 1. The fitting result was again acceptable. Also, the escaping photon 
effect has a relatively little effect on the fitting parameters except the normalizations. All 
the parameters are consistent with the previous ones before including this effect within 90% 
errors. The absorption-corrected X-ray luminosity is 2xl0 31 ergs s -1 . 

In addition, we also conducted the same spectral analysis to a C2 region, the soft X- 
ray clump seen in figure 1 (c). The obtained photon index and temperature are similar to 
those of the CI region within errors, except a lower absorption column density (0.1 x 10 22 
cm 2 when fitted with the leaky absorption model including the escaping photon effect). No 
sign of emission lines was seen. Hence, the spectral hardness of the extended emission is 
considered to be common within the whole region, and the offset peak in the soft band map 
is simply a consequence of a slight reduced absorption. 



4. Luminosity Function 

Based on the surface brightness of the extended emission, we estimated contribution 
from unresolved faint point sources to this emission, in order to know whether it can be 
explained by faint sources that are individually undetectable. We followed the same way of 
Ezoe et al. (2006) in which we utilized the luminosity function of the detected point sources. 
We below utilize the X-ray surface brightness of the extended emission after subtracting the 
escape photons from the point sources, obtained in the last section. 

Figure 4 shows the luminosity function of the 176 point sources in the CI. We derived 
the absorption-uncorrected 0.5-7 keV luminosity L. The X-ray flux of each point source is 
obtained by spectral fitting for a bright source (> 30 net counts or counts after subtracting 
the background), while by using a count-to-flux conversion factor derived from the summed 
point sources 1.6xl0~ n ergs s _1 cm s -2 (net counts/s)" 1 for a fainter source. 

The source number density increases toward lower luminosities and saturates below 10 
net counts, corresponding to the completeness limit of this observation. We here estimated 
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the limit from source number histograms as a function of the logarithm of net counts, and 
regarded the maximum of the histogram as the completeness limit. When utilizing the above 
conversion factor, 10 net counts correspond to 4xl0 28 ergs s -1 , in terms of the absorption- 
uncorrected 0.5-7 keV luminosity. This limit is one of the most lowest ones among those in 
the past Chandra observations of MSFRs. 

The solid line in figure 4 indicates the necessary point source number, in order to explain 
the extended emission by point sources. The point number is clearly short, which supports 
that the extended emission is truly diffuse. Furthermore, even when we extrapolate the 
luminosity function toward a lower limit of ergs s -1 , using a linear function fitted in 
log-log space from the completeness limit to 10 30 ergs s _1 , the estimated contribution of 
unresolved sources is at most ~ 10% of the extended emission. Based on these results, the 
extended emission of NGC 2024 can be considered as diffuse in nature. 



5. Discussion 

We have found the diffuse X-ray emission in NGC 2024. Because of the proximity of 
NGC 2024 and the long exposure time, this discovery itself is one of the strongest supports for 
the existence of the diffuse emission in MSFRs, among the previous Chandra results. Also, 
this result provides us with a new observational evidence that, even in a MSFR where only 
late O to early B stars exist, the diffuse X-ray emission can be observed if the sensitivity 
is enough high. The spectral analysis suggests that the diffuse emission is dominated by 
the thermal emission. At the same time, since the continuum (kT > 8 keV) of the diffuse 
emission is harder than typical spectra of young stars, a part of the emission may come from 
non-thermal origin. In NGC 6334, the soft-thermal and hard possibly non-thermal regions 
can be spatially distinguished (Ezoe et al. 2006). In NGC 2024, the soft and hard regions 
may be co-spatial and hence be observed as a mixture of both thermal and non-thermal 
emission. 

We then discuss whether the thermal or non-thermal interpretation is feasible in terms 
of energetics. First, if we assume that the whole diffuse emission is thermal, the total plasma 
energy U will be ~ 10 47 ?7 a5 ergs s _1 where i] is a filling factor, from the equation (3) in Ezoe 
et al. (2006). Using the X-ray luminosity L x and U, the plasma cooling time is estimated as 
U/L x ~ 10V' 5 y r > whi ch is far longer than the age of NGC 2024, from 0.3 (Meyer 1996) to 
several Myr (Comeron et al. 1996). Although the sound crossing time in a 10 keV plasma 
across the region of 0.5 pc in size is ~ 10 3 yr and hence short, the plasma with the estimated 
pressure of ~ 5 x 10 7 ??~ a5 K cm -3 , calculated from the equation (5) in Ezoe et al. (2006), 
may be confined by the dense HII region known to exist around the molecular cloud dark 
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lane (figure 1 c) (Subrahmanyan et al. 1997) and also the strong magnetic field within the 
molecular cloud (Crutcher et al. 1999). Then, by dividing the total energy by the assumed 
age of NGC 2024 of 0.3 Myr, the average energy input is estimated as 10 34 ergs s _1 

Second, we consider the non-thermal possibility. The flat continuum (r ~ 1 or kT > 8 
keV) of the emission strongly suggests the bremsstrahlung emission by 10 keV to several 
MeV electrons, rather than the synchrotron or inverse Compton emission (Ezoe et al. 2006). 
Since the Coulomb loss overwhelms the bremsstrahlung emission, if we assume that the 
diffuse emission is totally non-thermal, the necessary kinematic energy to be supplied is at 
least ~10 5 times larger than the observed X-ray luminosity, ~2xl0 36 ergs s _1 . 

One of the most plausible energy sources for the diffuse emission is the shocks generated 
by fast stellar winds (~ 2000 km s -1 ) from young massive stars. Its huge kinematic energy of 
the stellar winds can be easily converted via the strong shocks, among dense molecular clouds 
and HII regions, into the thermal (and non-thermal) energy of particles in the surrounding 
gases. This explanation has been proposed in M17 (Townsley et al. 2003) and NGC 6334 
(Ezoe et al. 2006). In NGC 2024, there are at least one late O to early B star IRS 2b (08V- 
B2V), one early B star candidate IRS 2, and 7 possibly massive protostars. As shown in figure 
1, these massive stars are spatially associated with the diffuse emission. A typical kinematic 
luminosity of the stellar wind is ~ i0 33-35 er g S g- 1 p er one B2V-08V star (Howarth et al. 
1989; Prinja 1990; Panagia 1973). Therefore, the necessary energy supply of the thermal 
interpretation is explained if there is at least one massive star with a strong wind, comparable 
or stronger than those of a typical B0.5V star. The non-thermal interpretation is also possible 
if we consider all the 9 sources are massive stars and have strong winds. Wind-wind collisions 
may effectively increase the energy of a shock. Hence, the stellar wind scenario is possible 
from the viewpoint of the energetics. 
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Fig. 1. — Adaptively smoothed X-ray images of NGC 2024 taken with the Chandra ACIS-I. 
Panels (a) and (b) correspond to 0.5-2 keV and 2-7 keV band images before applying the 
point source mask, respectively, while panels (c) and (d) are point-source excluded contours 
in 0.5-2 and 2-7 keV bands overlaid on the optical DSS image. Co-ordinates are the J2000. 
All the images are corrected for the exposure and vignetting, but the background is not 
subtracted. Boxes show the ACIS-I fields of view, while large circles (CI and C2) are 
regions utilized in our spectral analysis. Small circles are positions of IRS 2b, IRS2 and FIR 
1-7. Because IRS 2 is located just 5" south-east of IRS2b, these two circles are overlapped. 
The color intensity is plotted logarithmically from 5xl0 -10 to 5xl0 -6 counts s -1 pixel -1 
cm -2 in panel (a), while from 7xl0 -10 to 7xl0 -6 in panel (b). Also the contours are plotted 
logarithmically from 7.5xl0 -10 to l.OxlO -9 counts s -1 pixel -1 cm -2 in panel (c), while from 
1.3xl0 -9 to 4.0xl0 -9 in panel (d). 
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Fig. 2. — (a) ACIS spectrum of the CI region compared with that summed over 176 point 
sources, (b) The CI spectrum fitted with a power-law plus a narrow Gaussian model, (c) 
The same as panel b, but for a single temperature plasma model suffering two absorption 
column densities (solid and dotted lines), (d) The same as panel c, but including the escaped 
photons from the summed point sources (red lines). 
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Fig. 3. — The summed point-source spectrum in the CI region, fitted with the leaky absorp- 
tion model (solid and dotted lines). 
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log (L [erg/s]) 

Fig. 4. — Luminosity function of the point sources within the CI. The vertical axis is the 
column density of the source number, while the horizontal axis absorption-uncorrected 0.5-7 
keV luminosity. Errors are la Poisson. The thick line indicates the number of putative 
point sources at a given flux that can account for the observed total surface brightness 
of the extended emission. Three short vertical lines are the typical X-ray fluxes of point 
sources with 10, 20 and 30 net counts. The dashed line indicates the expected number of 
extragalactic sources, estimated from Giacconi et al. (2001). 
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Table 1: Parameters of the spectral fitting to the extended emission in the CI region. 



Model a 


N H i b 


A f H2 b 


T or kT c 




Norm.l° 


Norm.2 c 


S C Norm -finc 




Lk 


X Vd.o.f. 


po+ga 
leaky abs.l 


1 o + U.b 

0.62(< 1.1) 




11+8.8 
—3 5 


1 + I.6 


i +U.b 
i "^-0.8 
2.2±i"? 


9 2+ 2 ' 7 


u -'-0.2 u -°-4.0 


1.1 
1.1 


2.9 
4.0 


99.2/65 
85.9/64 


leaky abs.2 


0.21(< 1.1) 


3.3tl; 


11(>7.6) 


2.5(>1.0) 


0.73lg';S 


4 4+1-6 




0.63 


2.2 


84.2/64 



a Fitting models. The po+ga, leaky abs.l and leaky abs.2 indicate the power-law plus a 
narrow Gaussian model, the leaky absorption model without and with the escaping photon 
effect from the point sources, respectively. 

b Interstellar absorption, with N u being the hydrogen column density in 10 22 cm -2 . In the 
case of the leaky absorption model, two column densities are given. 

c T is the photon index in the power-law model, while kT is a plasma temperature in keV. 
d Metal abundance in solar, used in the leaky absorption model. 

e Normalization is photon flux at 1 keV in 10~ 4 photons cm~ 2 s _1 in the power-law model, 

and 10 _18 /47rD 2 EM in the leaky absorption model, where D is a distance to NGC 2024 

and EM is an emission measure in cm~ 3 . 

f A line center energy in keV. 

g Line intensity in 10~ 6 photons cm -2 s _1 . 

h The 0.5-7 keV flux in 10~ 12 ergs s" 1 cm s~ 2 . 

1 The absorption-corrected 0.5-7 keV luminosity in 10 31 ergs s _1 . 
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Table 2: Parameters of the spectral fitting to the summed point sources in the CI, derived 
by the leaky absorption model a . 



param leaky abs. model 



iVH2 

kT 

7 b 
7 b 



Zs h 

7 b 
^Ar 

7 b 

7 b 
^Fe 

Norm. 1 
Norm. 2 

Fx 

L x 
X 2 /d.o.f. 



0.55±0.09 



-0.3 
-0.2 



2.7" 

4.41°;! 

2.9(> 1.2) 
1.1±0.8 
0.06(< 0.32) 



0.93 



hO.37 
-0.36 



1 2 +0 - 7 

0.59(< 1.2) 
n 90+0.04 

u - zo -0.03 

17±« 
130±10 
9.5 
43 

292.1/157 



a Notations and symbols are the same as the leaky absorption model in table 1 except 
abundances. 

b Metal abundances of individual elements in solar, while those of other elements are fixed 
at 1 solar. 



